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The screening nature of the potential between external quarks in 
massless SU(N C ) QCD2 is derived using an expansion in Nf- the 
number of flavors. Applying the same method to the massive model, 
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1 Introduction 



Large N c expansion]!], and strong coupling expansion were shown to be 
powerful tools in analyzing various aspects of QGD^. In the present work 
we demonstrate the usefulness of yet another technique, that of expanding in 
large number of flavors Nf. We apply this method to derive further evidence 
for the screening behavior between external quarks in the massless theory, 
and the confining one in the massive theory 0. 

The screening nature of the potential was argued in ref.|| by substituting 
a static (abelian) solution of the equations of motion for the gauge configura- 
tion into the expression of the potential. This approach falls short in proving 
the screening nature of the interaction between the external quarks because 
of the following drawbacks, (i) Classical configurations cannot always recast 
for the quantum behavior of a system, (ii) Treating the external quarks as 
classical sources is strictly speaking justified only for large color representa- 
tions in a similar way that only particles with large angular momentum can 
be described in quantum mechanics as c-number quantities. 

In the present paper we suggest improvements that enable one to over- 
come these obstacles. To justify the fact that functional integral is dominated 
by the saddle-point configurations we introduce Nf flavor degrees of freedom 
and analyze the system in the limit of large Nf. It is easy to realize, espe- 
cially in the bosonization formulation, that the action can be brought into a 
form where it is multiplied by an overall factor of Nf. This obviously implies 
that j^- plays the role of % so that the limit of Nf — > oo corresponds to the 
classical limit of h — > 0. 

Moreover, the abelian gauge configuration that solves the equations of 
motion and is used in || to derive a screening potential, can be used as a first 
iteration in a systematic expansion of the gauge configuration in powers of 
One can thus check whether the screening nature persists also in corrections 
which are proportional to higher powers of 4-. However, if one presents 

the external quarks in terms of commuting functions as was used in [[|, 
then the iterated gauge field solution at any order will remain abelian if the 
leading one is. By abelian we mean that commutator terms in the equations 
of motions vanish. One may be suspicious that in that way one abelianize 
QCD 2 in an unjustified manner and hence the form of the potential cannot 
be faithfully examined. This situation can be avoided either by searching 
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genuine non-abelian solutions of the equations of motions [|[ or by using a 
different method of representing the external quarks so that the non-abelian 
commutator terms do not vanish. A proposal of this nature was introduced in 
refs. [0, where the density of the external quarks was represented in terms of 
non-commuting matrices. We compute the potential including corrections of 
next and next to next order corrections and show that it tends to a constant 
for large separation distance, thus it exhibits a screening behavior. 

Two dimensional massless QCD admits a universal nature 0. The intro- 
duction of mass terms to the quarks creates a much more dramatic alteration 
in the system than in the real world 4D QCD. This property manifests itself 
in the baryonic spectrum as well as in the form of the potential between 
external quarks. External quarks in the fundamental representations that 
were screened by massless adjoint dynamical fermions, find themselves in 
the confining phase once the dynamical fermions acquire mass. In ref . this 
phenomenon was derived for the analog abelian well as for QCD 2 

with two and three colors. In Q it was further shown for a general SU(N C ) 
gauge group by bosonizing both the dynamical and the external quarks and 
the conditions for finite energy quark configurations were analyzed. An ex- 
plicit determination of the potential, in a similar manner to way it was done 
in the massless case, is technically much more tedious due to the form of 
the bosonized mass term. Once again the expansion in large Nf comes to 
our rescue. The mass term can be expanded in terms of a (non-local) power 
series of the color current divided by Nf. Restricting the expansion to the 
lowest power of jj- renders the equations of motion in terms of the currents 
to tractable ones. Using this method we show that the leading solution for 
the potential exhibits a confining behavior. 

A natural framework that "embeds" the QCD 2 model with adjoint 
fermions is that of N — 1 super Yang-Mills theory. In ref. [§] the potential 
of that model is conjectured to be a screening one. We show that a naive 
bosonization of the fermionic degrees of freedom of the model leads to a 
wrong picture. It is further shown that by considering the contribution of 
the scalar field to the one loop vacuum polarization one indeed discover the 
non-confining potential. 

The paper is organized as follows. In section 2 the equations of motion 
of multi-flavor massless QCD 2 in the presence of external quarks are derived 
in the A_ = gauge. A large Nf iteration analysis is described in section 
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3. The potential is shown to be a screening one even for a non-commuting 
presentation of the external density. Section 4 is devoted to the massive 
multi-flavor case. Expressing the mass term in a power series in the currents 
we derive the confining potential for that model. An analysis of the potential 
for the N = 1 supersymmetric Yang-Mills theory is presented in section 5. 
We show that a one loop effect in the scalar sector which is the analog of 
the anomaly term for the fermions combines with the latter to produce a 
screening potential. 



2 Mult i- flavor QCD2 with external source 

Massless multi-flavor QCD2 with fermions in the fundamental representation 
of SU(N C ) and external current coupled to the gauge field is described by 
the following action 

S = J d 2 x tr{-^F% + if - eA,£ xt ) (1) 

where \I/ = i = 1 . . . N c and a = 1 . . . Nj and the trace is over both the 
color and flavor indices. Variation with respect to the Dirac field and gauge 
fields leads to the following equations of motion: 

= (2) 
D,f ext = (3) 
D,F^ = e(r + f ext ) (4) 

where = <9 M + e[A^, }. The first (0) and the second (||) equations are 
the covariant conservation of the dynamical and external vector currents, 
respectively. In addition the chiral anomaly equation reads 

eNf 

Dfl = ^T 6 ^ (5) 

where J 5fl = vj/^^vl/. The anomaly equation can be obtained by the equa- 
tions of motion of the bosonized action of (|l|) (see for instance ||). 

In order to maintain gauge invariance of the action (0), the external 
current should be covariantly conserved. This implies that, in contrast to 
the abelian case, one cannot set j® xt as a time independent function and 



3 



then fix j\ xt to be zero. Instead we fix the value of j^ xt and treat j\ xt as a 
"dynamical" variable of the problem. 

By choosing the A_ = gauge, using light cone coordinates x ± = -^(x°± 
x 1 ), and the two dimensional property J 5m = —e^Jy, we translate equations 
d),(D and © into 

d+J + + d-J- + e[A+, J + ] = (6) 

eNf 

d + J + - d_r + e[A+, J + ] = ^d^A + 
-d 2 A + = e{J++ 3 t xt ) 

In our analysis we will be interested in static solutions of the equations 
that correspond to static external sources. For this case the set of equations 
combined with the external source constraint takes the following simplified 
form 

^ dl J+ + e[A + ,n = -^l= dl A + (7) 
\8lA + = -e { J+ + ±= P +±=j) 
dxj + -^e[A + , p] + ^e[A+J] = 

where we have used the notation j = j\ xt and p = j® xt . Given p, equations 
(0) can be used to determine the dynamical variables A + , J + and j once 
boundary conditions are specified. 

The gauge field A + itself, is not the physical quantity of interest, however, 
once it is determined, the potential energy between the external charges can 
be found by substituting it in the effective action. Using the equations of 
motion the potential takes the following form: 

V = \ejdx tr{A +J t xt ) = \e J dxA\±={p + jf (8) 

where a are the adjoint indices of the color group a = 1, N% — 1. As was 
discussed in the introduction, the behavior of the potential as a function 
of the distance between the external quark and the anti-quark determines 
whether the theory is screening or confining. A linear potential means a 
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constant force, namely, confinement, while a constant potential (at large 
distances) means screening. 

We are now facing the question of how to incorporate the external sources. 
Consider a system of a quark in the fundamental representation placed at a 
distance of 2R from an anti-quark that transforms in the anti-fundamental 
representation. This can be expressed as the following classical c-number 
function 

p a = 8 al {5{x-R) -8{x + R)) (9) 

Strictly speaking, one is allowed to introduce classical charges and neglect 
quantum fluctuations only if the external charges transform in a large color 
representation. (This is an analog of the statement that only for quantities of 
large angular momentum quantum fluctuations are suppressed.) Moreover, 
by choosing ([]), there is an obvious "abelian" self-consistent solution^] of the 
equations ©for which all the dynamical quantities (A + , J + and j) points in 
the '1' direction and thus all the commutators vanish. This solution does not 
reflect the non-abelian nature of the theory. This situation is not avoided 
also when one uses an iterative expansion in large Nf, as is described in 
the next section. Hence, it seems that deriving conclusions from the abelian 
solution based on the use of (|[) is unjustified and may miss the true nature of 
the interaction. One way to overcome these obstacles is to search for "truly" 
non-abelian solutions of the equations. This approach was followed in ref. 
M. Here we proceed by implementing Adler's0 semi-classical approach for 
introducing static external quark charges. In this approach the quarks color 
charges satisfy non-abelian SU(N C ) color algebra so that the external quark 
charge density takes the form 

p a = Q a S(x — R) + Q a 5(x + R) (10) 

Q a and Q a are in (N c , 1) and (1, N c ) representations of SU(N C ) <S> SU(N C ) 
group respectively. The algebra of those operators which was worked out in 
[| is reviewed briefly in the appendix. Note that unlike what followed 
from the classical expression (|9]), now since Q a and Q a do not commute, 
"non-abelian" solutions are expected. In the following section we apply a 
systematic 4- expansion of A + and derive such a solution. 
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3 Large Nf expansion and non-Abelian 
solutions 

The form of the set of equations (|7|) suggests a natural expansion. Large 

Nf, with ^— ^ = /i 2 fixed, means weak coupling constant e. In this case a 
procedure of extracting a solution by iterations can be applied in the following 
way. A solution for A + , J + and j of the equations expanded to a given order 
in e is inserted back to (|7|) as a source to determine the next order solution. 
A similar treatment in four dimensions is given in refs.|8], |j. 
The formal expansion in e is as follows 

A + = eA^ + e 3 A^ + e 5 A^ + ... 

J+ = j(0) +e 2 J (2) +e 4 J (4) + ___ (n) 

J= e 2 j( 2 )+e 4 j( 4 ) + e 6 j( 6 ) + ... 

Substituting A + , J + and j back in (|7|) one gets equations for A^\ J^" 1 ) and 
which to lowest order in take the form 
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N 



f 



J_ dlJ (o) = _ l _^J_ dA ii) 
\dlA^ = -(JO) + -Lp) 
Assuming vanishing currents at infinity the equations can be rewritten as 

J(0) = _I M 2 A (1) 



(«9 2 - p 2 )A« = -V2p 



The solution for A^ is 



^(i) = |5(Q e -/^-*l + Q e -^+R\) (12) 

This is the "abelian" solution of Q] with a replacement of the c-number 
charges with the non-commuting ones. Whereas in the case of QED2 coupled 
to external sources, this is (with c-number charges) an exact solution, in the 
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present case it is the leading (in -^-) contribution to A + . The next to leading 
order set of equations is 

^d 1 J^ + [A^\j^] = - 1 -^d 1 A^ 

1%A® = -(J® + ±j<*) 
d l3 ^ + ^=[A^\p)=0 
Substituting A^ and we find 

J(2) = -1^(3) 

(d 2 - p?)A^ = -yfij® 
with j( 2 ), which is determined by the previous iteration, of the form 

J {2) = P] = ^[Q, Q]e-^ R (e(x + R)- e(x - R)) (13) 

where e(x) is the step function e(x) = 1 for x > and e(x) = —1 otherwise 
and ^j- denotes the integral Jf ^ dx' . 

The latter expression acts as a source to A^ which leads 



A^ = ^[Q,Q}e- 2 » R {e(x + R){l-e-^ x+Rl )-e{x-R){l-e^ x - Rl )) (14) 



V2, 

Far from the sources A*- 3 -* approaches zero. This indicates that the poten- 
tial is approaching a constant value when one quark is taken to be far from 
the other. However, as we shall see, A^ will be needed to observe the first 
correction to the potential. This calculation is written in the Appendix. 

We substitute now the expression found for A + and j + in (||) to derive 
the following potential 



If 1 
V = -ej dxA + -j=(p + j 



J dx(eAM + e 3 A^ + e 5 A^ +...){p + e 2 / 2) + e 4 j {4) + . . .) 
^ f dx (e 2 A^p + e 4 (A (1 ¥ 2) + ^ (3) P) + e%A^j^ + A^j® + A®p) + . . .) 



2^2 
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All the e 4 terms vanish because they contain [Q,Q}Q ~ f abc Q a Q b Q c part 
which vanishes. Therefore the first non-trivial correction is 0(e 6 ). 



e 2 



V(2R) = — {QQ + QQ + 2QQe- 2 » R ) (15) 



[Q, Q] 2 \(l - e~ 2flR ) 2 - 2fiR e~ 2 ^(l - e~ 2 ^) 



8/i 5 ' 

Let us compute the group constants that appear in (fL5|) . QQ stands for 
Q a Q a which is the second Casimir operator of the fundamental 
representation- Q a Q a = . Similarly QQ = . The value of QQ is 
determined easily by using 2Q a Q a = {Q a + Q a ){Q a + Q a ) - Q a Q a - Q a Q a . 
In the singlet coupled state Q a + Q a = and hence Q a Q a = — 2 C JV T 1 . 

— 2 — 

The value of [Q, Q] is computed by the QQ algebra (see the Appendix), 
itis[Q,Q] 2 = f^-i. 

Defining d = 2R , the potential takes the form: 

+u( ) f (1 - e-^y - fid e-^Cl - 6-"*) 

^ K 2N f J 2 2N C " J P v ; 

Thus, the potential that includes the first correction to the abelian one 
approaches a constant value at large distances where the force between the 
external quark and the anti-quark vanishes. The obvious question now is 
whether we can infer from this result that the screening nature of the inter- 
action remains valid to all orders in We cannot provide a general proof 
of that statement, but we believe that indeed that is the exact nature of the 
interaction. This is based on the following argument. The structure of all 
higher contributions is (df — [i 2 )A = —j, where j is determined by previ- 
ous iterations. Therefore, j would always vanish far from the sources and 
consequently the gauge field A would exhibit the same behavior. 



4 Large Nf expansion of massive QCD2 

Let us consider now the case of massive dynamical quarks. Whereas, for 
the massless case the equations that determine the potential can be derived 
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using both the fermionic picture and the bosonized one, here we can apply 
our analysis only in the bosonization description. The bosonized action of 
massive QCD 2 with Nf fundamental representations in the gauge A_ = 
takes the following formpj 



d x tr 



m 2 : (ghlexpi-i^j^-fy + exp(^^^-0)/ t /i t ^ t ) + 



\{d-A + f - eA + (J+ + j+ t ) 



47T 

N c N f 



where tr is over U(N f x N c ), g G SU(N f ), h G SU(N C ), exp(- 
17(1), I G U(N C x N f )/SU(N C ) x SU(Nf) x U(l), S$f w is the level k WZW 
action. : :^ denoted normal ordering at mass scale m and m 2 = m q mC 
where m q is the quark mass, and C = |e 7 with 7 Euler's constant. 

The bosonized mass term is the only term in the action that couples the 
colored and flavored sectors. Since we are interested only in the form of the 
potential we can simplify the analysis of the equations of motion of the full 
theory by restricting ourself only to the colored sector. This can be achieved 
by setting g — 1, I = 1 and = 0. The normal ordering can be performed 

first at the scale fj, = 7^ and then one can replace the two mass scales 
H and m q by a single scale by normal ordering at a certain scale m^. In 
that case m = [N fm q C '/i 1_Ac ] 2 - A c where A c the dimension of h is given by 

^ c N C (N C +N f ) 

which leads the following action 



S = S^ w (h) + J d 2 x tr \m 2 (h + h- 1 ) + \{d-A + f - eA + (J+ + 



(17) 

and the trace is over the color degrees of freedom. Equations (|6]) now 



read 



eN 



d + J + + e[A+, J + ] = — } -d^A + - im 2 {h- 1 - h) 



2tt 



-d 2 A + = e{J + +j 



ext ) 



(18) 
(19) 
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The difference from the equations derived in the massless case is that 
here there is a dependence on h on top of the dependence on the currents. 
Unlike the abelian theory here we can write a closed expression for the group 
element (h) only in a formal form in terms of the current. However we can 
express it as an expansion in powers of J + . Recall tht the current J + is 
related to the color group element h by 

j+ = l ll hd _ h -^ (20) 
2n 

The expanding of the inverse relation takes the form 

ft " 1 = 1 + ^ + + (^' 2 ^ + ^ + ) + '" < 22 > 

where the dots denote additional terms which are higher powers of J + . This 
expansion makes sense upon the substitution of J+ = £^° =0 )"J (2n) given 
in eqn. ([TT|) . Moreover, note that even for the free level Nf WZW model this 
expansion is justified since J + behaves like J + ~ \j~N~f as can be deduced 
from the associated affine Lie algebra. For the expression needed in (|T8| ) we 

Thus, the set of equations (0), for static solutions, take the following form in 
the presence of mass 

_ !rf( _ 2 ^i J++( ^) 2 i|i J+ , J+1) 

v iN f d x K iN f 1 d x [ di ' u 
ld 2 1 A + = -e(J + + -^p+-^j) 

dij + -jf\ A +ip\ + 7f e [ A +'^ = 
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which to leading order, in is 

1 n + eJV f 1 „ , 9 \/22tt 1 T , 

^ + = -e( J + + i=p) 

Eliminating J + , we find the following equation for A + 

e 2 N 2 

(! + TT^) 5 ?^ - - e ^P ( 24 ) 
The solution of the equation in the presence of ([H]) is 

e e 2 iV? 

A + = —={l + ^4)- 1 (Q \ x-R\+Q\x + R\) (25) 
Substituting A + in the potential yields, 



e 2 e 2 N 2 f 

V = --(l + -^-)~ l QQx2R (26) 



Using the definition d = 2i?, substituting m 2 , /z 2 and QQ which is minus the 
quadratic casimir operator of the test charges in a general representation 1Z 
we obtain 

^ = f^(l + ^-^)- 1 ^)xd (27) 
2Nf S7rCm q 

The same expression for the potential in the abelian case was obtained in 
[Q]. Thus the dominant contribution exhibits a confinement behavior. 
It should be emphasized that in the above analysis it is assumed that the 
external charges cannot be composed by the dynamical ones. This is the 
analog of the abelian case where the external charges are not integer multiple 



of the dynamical charges [10]. It is well known that in systems, where the 
latter does not hold, the string between the external quark and antiquark 
can be torn apart by a pair creation. Similarly in the non-abelian case we 
expect that when the test charges can be composed of a multiplication of the 
dynamical charges the string could be torn. Confinement is restored (and 



(27) is valid) only when when the test charges cannot be composed of the 
dynamical charges |TT|. It is clear that in the latter case the confining nature 
of the theory survives higher order corrections. 
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5 Supersymmetric Yang-Mills 



Next, we would like to investigate the large distance behavior of N = 1 two 
dimensional supersymmetric Yang-Mills theory, and test the conjecture that 
the system is in a screening phase ||. 
The SYM2 action [|12|] is the following 



S = Jd 2 x tr + ^\$>\ + \{D^) 2 + 2ief\ lb \} , (28) 

where the gluon field, A the gluino and <fi a pseudo-scalar - the compo- 
nents of the vector supermultiplet- transform in the adjoint representation 
of SU(N C ). 

The action is invariant under the following supersymmetric transforma- 
tions 



8A ll = ie7 5 7 M v / 2A 
5(f) = -eV2X 

% 



where e is the fermionic parameter of transformation. In order to determine 
the behavior of the theory it seems natural to follow the same procedure 
used in the previous sections, namely, to bosonize the fermionic degrees of 
freedom. Then, to proceed by solving the equations of motion for the gauge 
fields and deducing the potential between the external sources. In the present 
case one needs to bosonize the gluinos which transform in the SU(N C ) adjoint 
representation. The bosonized version of adjoint fermions can be expressed 
in terms of a SU(N C ) WZW model of level iV c [0]. The form of the Yukawa 
term follows from the identification of : h — h^ 1 : with : A75A :. The full 
bosonized action in the A_ = gauge takes the following form 

Q qWZW (OOl\ 

^bosonized — D (N C ) l zy J 



+ J d 2 x tr (\(d~A + ) 2 + ^(D^) 2 + ifhecf) : (h - h' 1 ) ^ -eA + (J 



1 Jext, 



where rh is the scale at which the normal ordering is performed in a sim- 
ilar way to the one introduced in eqn. (|17D . Eliminating the scalar and 
bosonized fermion degrees of freedom from the equations of motion which 
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arise from(p9|), one finds that the effective equation of motion of static gauge 
field is 

The meaning of this equation is that the corresponding propagator of the 
gluon has three poles and one of them is tachyonic. It turns out, as will 
shown below, that the source of this unwanted pole is the fact that the 
fermionic contribution includes quantum corrections in the form of the axial 
anomaly, whereas the scalar contribution is a tree level one. This follows 
from the fact that we write down the classical equations of motion of (^) 
but it is well known that the tree level bosonized action incorporates the 
axial anomaly. 

The correct procedure for the supsrsymmeric case is to consider loop 
effects of both the fermions and the scalars on the gluon propagator. 

Before we describe the calculation associated with SYM 2 model, let us 
review some properties of the Schwinger model and two dimensional scalar 
electrodynamics. 

The lagrangian of massless QED 2 (the Schwinger model) is 

C = ~F% + M0* (31) 

The resulting classical equation of motion for the gauge field is 

DF = ee» v d„j v , (32) 

where F = -e^F^ and j M = ^7 M \1/. Classically, the right hand side of the 
above equation vanishes because d^e^ j u = —d^jg = 0. However, quantum 
mechanically the axial current is anomalous d^jt^ = -F and therefore the 
quantum mechanical form of equation (|32|) is 

e 2 

(□ + — )F = (33) 

7T 

which leads to the screening behavior. 

In scalar electrodynamics similar phenomenon occurs. The lagrangian of 
the theory 

£ = -^ + (^)W (34) 
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yields the following classical equation of motion for the gauge field 

DF = ee^d, 3u , (35) 

where the scalar abelian current is = —i (0*<9^0 — 0(9^0* — 2ieA M 0*0). 
Note that unlike in the fermionic sector, the right hand side of (|35"D does not 
vanish even at the classical level since there is no conserved axial current in 
the scalar sector. The classical divergence of the "axial current" is 

d^j v = -ie^(2d^d u <j) - 2ie<9 M (A,0*0)) (36) 

This relation is modified quantum mechanically similarly to the modification 
in the Schwinger model. The one loop vacuum polarization diagram modifies 
the axial current non conservation. In scalar electrodynamics the one loop 
vacuum polarization is given by the following two Feynman diagrams (the 
dashed line is the scalar field) which combine to 

IMP) = (^ - ^) x e - f dx , 2 ^ (3 \ f) 2 (37) 
k 7i Jo k z x{± — x) — m z 

This expression suffers from severe infra-red divergences in the m 2 — > 
limit IP- 4||. This technical obstacle can be overcome by the use of a point 



splitting of the vector current [15 



f = lim exp(-ie / A^dy^x 

e-^U J x — ^ 

-i (>(x + e -)dy(x - e -)- </>{z - 1 W(x + e -)- 

ieA»(x + ^)</>(x - € -)^{x + ~) - ieA»(x - |)0(x - ^)0*(x + ~ 

In fact, the point splitting method enables us not only to get rid of the infra- 
red divergence but also to derive the quantum mechanical modified version 
of© 

d^ju = ~F - ie^{2d^d v( f> - 2zefy(4,0») (38) 
lix 



which leads to 



e 2 



'a + —)F = -ietr&dnPW - 2ied ll (A v <f>*<f>)) (39) 

Z7T 
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Although the right hand side of the above equation does not vanish, it is 
clear that the nonlinear terms are interaction terms that can only modify the 
photon mass but can not make it massless. It is thus clear that the photon 
of scalar electrodynamics behaves in a similar way to the one of Qi?Z?2[0- 
Returning back to the SYM2 case two additional complications are in- 
troduced, (i) The gauge interaction is a nonabelian one and (ii) The gluino 
interacts with the scalars via a Yukawa term. Nevertheless, we can follow 
the same procedure used for the abelian models and compute the mass of 
the gluon pole. The gluon equation of motion now reads 

D^F = e^D^Jt + 4), (40) 

where denotes the gluino vector current and denotes the scalar vector 
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current. The equation for the divergence of the fermionic axial current. 

eN 

e^DJl = -F + ie(0AA + AA0) (41) 

7T 

The factor N c which appear in the anomaly term is due to the adjoint gluinos 
which run in the anomaly loop. A similar equation holds for the scalar current 

eN 

e^D.Jt = — —^F - e^(-2i[0, d v <f>] + 2<#, [A v , 4>]]) (42) 
Ztt 

Thus the quantum version of equation ( pE0|) is 

3p 2 7V 

+ -%f)F = ie 2 {f\\ + AA0) - ee^(-2i[0, d v (j>] + 2e[tf>, [A u , 0]]) 

(43) 

which means that the gluon propagator has only a single pole which is mas- 
sive. Though our results are based on one loop calculations, it seems that 
higher order corrections - which involve gluino and scalar interactions can- 
not spoil the massive nature of the gluon (but may shift its mass). The 
implication of the last equation on the potential between external charges is 
clear. The interaction mediated by the exchange of these massive modes is 
necessarily a screening one. The potential takes the form of eqn.flTED with a 
range that behaves like ~ [e\fN\\~ l . 



6 Summary 

The study of the interplay between screening and confining phases is a basic 
question in strongly interacting systems. Usually due to the lack of adequate 
perturbation expansion, exact statements about confinement versus screen- 
ing cannot be made. One system that evades this fate is QCD2. Special 
powerful techniques that are applicable only in 2D systems enable one to 
derive significant results. It is because of these results that QCD 2 may serve 
as an important laboratory to study real life QCD. 

Evidence that massless QCD 2 , regardless of the quarks representations, 
is in a screening phase was presented in the paper of Gross et. al.||. The 
dynamical quarks were shown to screen external charges even if the latter 
are in a representation that cannot be composed of those of the dynamical 
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ones. It was further shown J| that once a mass term is turned on the dy- 
namical fermions develop a non-vanishing string tension namely, a confining 
potential. Tensionless strings occur only in the particular cases that the rep- 
resentation of the external quarks could be gotten by a composition of the 
dynamical ones. These results were argued using several different methods. 
In particular, the potential was extracted by substituting an abelian solution 
of the equations of motion of both the abelian theory and certain non-abelian 
ones. 

In the present paper we derive additional supporting evidence for the 
picture drawn in ||: (i) We justify the use of the equations of motion by ap- 
plying a large Nf expansion; (ii) We extend the prove of screening mechanism 
by using "semi-classical" external charges; (iii) Show the confining potential 
for a bosonized massive model; (iv) Show that the 2D super YM theory has 
a screening behavior. 

A given massless multi-flavor QCD 2 model is a point in the two dimen- 
sional (N c , Nf) grid. The domain of Nf = 1 and large N c was described in 
the seminal work of 't Hooft in the form of the confining mesonic spec- 
trum. The analysis of 't Hooft was insensitive to the question of whether 
the quarks are massless or massive. One may wrongly get the impression 
that the massless theory confines. In fact, this limit is not adequate for 
the study of the question of confinement versus screening. Assume that the 
potential is of the form of leading term of fllpp. Recall that in the large 
N c limit, e 2 N c is kept finite. This implies that the potential behaves like 
(1 _ e ~vk R ) ^ -fe.R x (1 - \^f=R + o(^) ) for fixed R and large N c , with 

finite constant. Now it is clear that in the limit of N c — > oo 
the potential looks like a linear potential, and thus one cannot discriminate 
between the two scenarios. 

The opposite corner in the grid of theories is that of finite N c and large 
Nf with e 2 Nf kept finite. It can be shown [I7J that this limit corresponds to 
an approximate system of N 2 — 1 abelian theories. In that case the screening 
nature can be attributed to an exchange of Schwinger-like massive modes. 
The present work, as well as @, indicates that no "phase transition" should 
be expected when passing to models with small number of flavors. The 
screening nature of the potential may indicate that the theory with finite N c 
has in its spectrum states of masses of the order of e. In the large Nf, using 
the analogy with the massive Schwinger model, one can get a general picture 
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of the passage to a confining behavior. The mass of the massive state of the 
Schwinger model is shifted once quark mass is turned on. But an additional 
light state emerges. Exchange of the latter mode causes confinement. 

We would like to emphasize again, that studying the quantum system by 
analyzing the corresponding equations of motion is a justified approximation 
only provided that the classical configurations dominate the functional inte- 
gral. This condition is obeyed in the large Nf limit since in that case an Nf 
factor that can be put in front of the whole action of the colored sector, and 
thus plays the role of 

Another improvement over the analysis of is achieved by implementing 
the idea of || || to introduce the density of the external quarks in terms of 
non-commuting matrices. In that way the non-abelian nature of the large 
Nf limit of theory manifests itself in the form of non vanishing commutator 
terms whereas using the ansatz of || for the external charges has an abelian 
behavior in the same limit. 

It may seem that the analysis of the question of screening versus confine- 
ment in SYM2 could follow very similar lines as those of QCD2. However, it 
turns out that the equations of motion that follow from the action expressed 
in terms of bosonized gaugino fields, lead to unacceptable conclusion. In 
particular it reveals a tachyonic pole to the gauge field. This situation oc- 
curred due to an unbalanced treatment of the gaugino and scalar degrees 
of freedom. It is only after including quantum correction also to the scalar 
fields, in the form of point split currents, that a meaningful result could have 
been extracted. The latter corresponds indeed to a screening phase as was 
conjectured in |||. 

Many open questions associated with the interplay between confinement 
and screening are still unresolved. One is the derivation of the string tension 
for massive dynamical quarks in any representation of color group and for 
any representation of the external quarks. It is speculated that this string 
tension is ~ fim, where m is the dynamical quark mass and [i is proportional 
to the coupling constant, whenever the test charges cannot be composed of 
the dynamical charges, whereas when the test charges can be composed of 
the dynamical charges - the string tension is expected to be vanish. The 
super symmetric YM 2 has further extensions. In particular the large N c limit 
of the N = 8 supersymmetric case has an important significance in relation to 
the matrix model representation of M theory. Last but not least is obviously 
the implications to 4D QCD. In particular an interesting question is whether 
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the 2D large Nf expansion sheds any light on the 4D one. 
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A Appendix - The QQ Algebra 

The Appendix is based on refs. || |6| . 

Regarding the components of the Gluon fields as operators, ordering am- 
biguities may appear in products such as / diA+A+. These ambiguities 
are resolved by the transcription 

[u,v} a = ^f abc (u b v c + v c u b ) (44) 

where a, b, c, . . . = 1, 2, 3, N 2 — 1 are SU (N c ) indices. The algebra is spanned 
by four operators: Q a ® 1, 1 ® Q a , f abc Q b Q c and d abc Q b ® Q c . Where 
Q a , Q a are SU(N C ) generators. 

Using the following identifications 

e a l = ^{Q a ®l + l®Q a ), (45) 

4 = 4t(Q a ® i - 1 ® Q a ) - 4 rfabc ^ ® Q c > ( 46 ) 

e^-^rVW, (47) 

ei = ( N ] ~ 4) (g a ® i - 1 (g g a ) + c/ abc g b ® g c (48) 



4iV 



It was shown that 



[e*, e,-]" = e^e£ i,j,fc = 1,2,3 (49) 
h,e 4 ] a = (50) 



and the set is orthogonal 



tr e«e b = (51) 



This means that the SU(N C ) algebra is actually reduced to a SU(2)®U(1) 
problem. 

Using the above algebra, it is easy to see that tr [Q,Q]Q = tr [Q,Q]Q = 
and tr [Q,Q][Q,Q] = §*r e%e% = f^=±. 
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B Appendix - Calculation of 

This Appendix contains the calculation of and which are needed to 
determine the potential in the case of massless dynamical quarks (Section 3). 
The next to next to leading order set of equations which determines 

is 

4=9i^ (4) + [^ {3 \ J (0) ] + J m ] = -l^^dU^ 
V 2 2 y2 

5u (4) + ^[^,p] + i=[A« ) / 2 )] = o 

After substitution of A (1 \ A^ and J (0) , J (2) we find 



J(4) = _ I 2 A (5) 

2 P 

(9 2 - /i 2 )^ 5 ) = -v/2j 



(4) 



which is very similar to the leading and next to leading order sets of equations. 
The source can be calculated by using the values of A^ and A^ 

^=-U lAi3) - p] -7ik lAmj2)] (52) 

the result of the calculation is 



,(4) = 



8/i 3 

1 - e -^ R (l-e- 2 » R )e(x + R)x[[Q,Q},Q} 



8/i 3 

1_ ( e ( x _ R ) _ e ( x + _ e -Mk-«l) + (e(a . + i2) + i)(i _ e - 2 ^)) x [[0, Q], Q] 



^ (e(x + i?) - e(x - i?))(l - e-^l^l) + (e(x - R) + 1)(1 - e" 2 "*)) x [[Q, Q], Q] 



8/i 
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The first two lines in the expression of arise from the commutator of 
and p, the other lines are due to the commutator of and j^. Using 
a tedious calculation yields the following expression of A^ 



V2 



8p 5 
V2 
~8p' 

V2 
'8p 5 

V2 
8p 5 



e -W(\ _ e - 2 » R )e{x -R)(l- e-^ x - R \) x HQ, Q],Q] 
5 e-^ R (l-e-^ R )e(x + R)(l-e-^ x+R \) x [[Q,Q],Q] 



x 



-2pRe^ x ~ R ^ 



x < -R 
-R < x < R 



-2(1 - e ^ R ) + 1(1 - e -^) e -M*-«) x > R 



x [[Q,Q],Q] 



±(1 — g-4^R^ e At(a:+_R) 



X 



-2 + (/ix + pR + |) e -M*+«) -R<x<R 

_^l e -inR e n(x+R) 

-2(1 - e^ R ) + 2pRe~^ x+I V x > R 



x [[Q,Q],Q] 
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